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ABSTRACT 

The quasi-stellar object (QSO)/merger Mrk 231 is arguably the nearest and best laboratory for 
studying QSO feedback. It hosts several outflows, including broad-line winds, radio jets, and a poorly- 
understood kpc scale outflow. In this Letter, we present integral field spectroscopy from the Gemini 
telescope that represents the first unambiguous detection of a wide-angle, kpc scale outflow from a 
powerful QSO. Using neutral gas absorption, we show that the nuclear region hosts an outflow with 
blueshifted velocities reaching 1100 km s~^, extending 2 — 3 kpc from the nucleus in all directions in 
the plane of the sky. A radio jet impacts the outflow north of the nucleus, accelerating it to even higher 
velocities (up to 1400 km s~^). Finally, 3.5 kpc south of the nucleus, star formation is simultaneously 
powering an outflow that reaches more modest velocities of only 570 km s~^. Blueshifted ionized gas 
is also detected around the nucleus at lower velocities and smaller scales. The mass and energy flux 
from the outflow are ^2.5 times the star formation rate and >0.7% of the active galactic nucleus 
luminosity, consistent with negative feedback models of QSOs. 

Subject headings: galaxies: evolution — galaxies: ISM — galaxies: jets — galaxies: kinematics and 
dynamics — quasars: individual (Mrk 231) 



1. INTRODUCTION 



Theory predicts that negative feedback in active galac 
tic nuclei (AGNs) regulates the growth of supermassive 
black holes and their accompanying starbursts Mat- 
teo et al ||2QQ5] [Narayanan et al.|[2QQ8| [Hopkins <^ Elvis 



201 op . On the scales of galaxies, radio jets power neutral 



and ionized outflows with significant mass and energy 
flux fMorganti et al.'"2005: "2007'; 'Holt et al.[ |2Q08( |Fu fc] 



l^ockton 2009; Holt et al. 2011). Wide-angle outFows on 
sub-kiloparsec scales are also c ommon i n radio-quiet sys- 
tems (e.g., [Crenshaw & Kraemer 2005). Precisely how 



and to what degree these phenomena are the negative 
feedback predicted by theory is an area of active study. 

The ultraluminous infrared galaxy (ULIRG) Mrk 231, 
the nearest infrared-luminous quasi-stellar object (QSO), 
is a unique system because it hosts several types of AGN 
outflows. It also combines many of the phenomena of in- 
terest in model s of galaxy mergers that then form power - 
ful QSOs (e.g., |Sanders et al.| 1988] [Hopkins et al.j2005[ ). 
The nuclei in Ivirk 231 have coalesced, but pro minent 



featur es of an interaction are still present (e.g., Surace 



et al. 1998|. It contains a power ful, unobscured AGN 



alongside an obscured starburst ( |Veilleux et al. 2009). 
Finally, several type s of outflows are present: pc and 



kpc-scale radio jets fCarilli et al. 1998 [Ulvestad et al. 
|l999); a low ionization, broad absorption line o utflow 
with blueshifte d velo cities up to 8000 km s ~^ (e.g., Bok- 
[senberg et al. 1977 Boroson et al. 1991), presumably 
arismg at pc scales; and a kpc scale outflow with veloci- 
ties exceeding 1000 km s~^. 

The properties of this third outflow concern us here. 
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Blueshif ted emission lines led to speculation of its ex - 
istence (Hamilton & Keel 1987 Krabbe et al. 1997), 
but emission lines are subject to mismterpretation, as 
blueshifted lin es could indicate inflow or outflow (e.g., 
Veilleux|1991 ). Using long slit observations, 'Rupke e t al. 
([2005a[) showed conclusive evidence that a >1000 km 
neutral outflow was present at galactocentric radii up 
to 3 kpc. More recently, this outflow was detected in 
the molecular phase at very similar velocities, with a 
marginal detec tion o f extended emission ( [Fischer et al. 



2010; Feruglio et al. 



2010). 



The structure and power source of this outflow have 
remained elusive. Though the high velocities suggested 
the influence of the A GN was helping to power the wind. 



Rupke et al. (2005a) were unable to determine whether 
the wind was powered by the jet or more widely-directed 
AGN energy, and what role the starburst played_. The 



high molecular gas veloci ties observed (Fischer et al. 
2010[ [Feruglio et al.[[2010[ ) also suggest AGN influence, 
but the spatial resolution of these observations made the 
power source equally difficult to determine. 

In this Letter, we present integral field spectroscopy of 
the outfiow in Mrk 231, probing both the neutral and 
ionized phases using Na I D 5890, 5896 A interstellar 
absorption and the Ha / [N II] emission-line complex. 
Using sub-kpc spatial resolution, we untangle the struc- 
ture and dynamics of the outfiow and show that jets, 
wide-angle AGN infiuence, and the starburst all play a 
role in moving large amounts of gas out of the central 
regions of this galaxy. 



2. OBSERVATIONS, REDUCTION, AND ANALYSIS 
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Mrk 231 was observed on 2007 July 7 UT with the Inte- 
gral Fiel d Unit in the Gemini Mult i- Object Spectrograph 
(CMOS; lAllington-Smith et al.||2QQ2l |Hook et al.|^QQ4 ) 

on the Gemini INorth telescope. (Jloud cover was mmi- 
mal, and image quality was very high (0^5 - 0^6). We 
obtained five exposures of 900 s each, at a position angle 
of 35° to ensure that diffraction spikes from the nuclear 
point source were not oriented along the north-south jet 
axis. Our dither pattern centered on the nucleus, with 
exposures alternating ±V!5 along one of the field axes. 
We used the integral field unit in two-slit mode and the 
B600 grating, yielding wavelength coverage from 5600 to 
6950 A and a spectral resolution of 1.80A at 6300 A. 

We reduced the data using the IRAF data reduction 
package provided by the Gemini Observatory, supple- 
mented with custom IDL routines. The 0^3 spaxels in 
the final data cube maximize signal-to-noise ratio while 
adequately sampling the seeing disk. The resulting cube 
contains 525 spectra and measures 6'.'3 x 7^5, with the 
long axis oriented northeast to southwest. 

We modeled the continuum and emission lines in each 
spectrum using UHSPECFIT, a suite of IDL routines 
that fits a contin uum and emission lines to spectra 
( Rupke et al.||2010 ). We fit the continuum at each point 
as a linear combination of the nuclear spectrum and a 
smooth host galaxy continuum. Where possible, we fit 
two Gaussian velocity components to the emission lines. 
One component represents the narrow, rotating veloc- 
ity component, and the second represents an underlying 
broad, blueshifted component. 

We model ed the Na I D a bsorption lines using the 
method of (Rupke et al. (2005b). One velocity component 
provided good tits throughout the cube. The nearby He I 
5876 A emission line was parameterized using the fit to 
other lines, but in almost all spaxels this line was absent 
or very weak. 

Figure [T] shows six example fits to the emission and 
absorption lines. They illustrate the high quality of the 
fits. 

Our line profile modeling is based on Gaussian veloc- 
ity profiles. As such, we define outflow velocities in this 
Letter based on the properties of the normal distribu- 
tion. We define negative velocities to be blueshifted and 
outflowing. 

'^50% = center of Gaussian profile 

(50% of gas has lower outflow velocity) (1) 

^84% ^50% - ^ 

(84% of gas has lower outflow velocity) (2) 

'^98% ^50% - 2cr 

(98% of gas has lower outflow velocity) (3) 

3. RESULTS 

3.1. Host Galaxy 

Our data reproduce known properties of the host 
galaxy in Mrk 231. Figure [2] shows the entire galaxy, 
with the GMOS field of view superimposed. Ha emis- 
sion traces young star formation, including the edg e of a 
prominent b lue arc ~5 kpc south of the nucleus ( ^Surace 
et al. 1998). The [N II] /Ha map reveals high excita- 
tion outside of star forming regions; high extranuclear 
excitation is common in ULIRGs ( ^Veilleux et al.,^1995^ 



Monreal-Ibero et alJ[2006). 

ih( 



I'he rotation of the central gas disk of Mrk 231 has 
been mode led using CO observati ons fBryant & Scov- 



1^ 1996; D ownes fc Solomon|1998 ), yielding a projected 



velocity amplitude of 70 km s~^ along the 90° line of 
nodes and a disk inclination of 10° — 20°. T his molecu- 
lar gas is conc entrated in a r ~ 1 kpc disk (D ownes fc 



Solomon'1998'). From CO and H I observations"! Downes 
STSolomon 1998{ |Carilli et al.|1998 ), we adopt a systemic 
velocity of z = 0.0422. This yields a spatial scale of 0.867 
kpc arcsec"^. 

The rotation of the ionized gas lines up with that of the 
molecular component, despite the difference in extinction 
of the two components (Figure Isl). Our data trace this 
rotation to larger radii than the CD observations. Devia- 
tions from the typical galaxy rotation curve are evident, 
but discussion of these features is outside the scope of 
this Letter. 

3.2. Outflow 

The second ionized gas component is blueshifted and 
much broader than the rotating component. Within 
1.5 kpc of the nucleus, the center velocity of the second 
component averages ('^^50%) = —150 km s~^, with larger 
velocities closer to the nucleus. Because these compo- 
nents are also very broad, with (FWHM) = 700 km s~^, 
the ionized gas reaches ('^98%) = —760 km s~^. The re- 
gion of broad, blueshifted emission is asymmetric, with 
higher velocities to the north and more extended emis- 
sion to the east. 

The neutral atomic gas in Mrk 231, as traced by 
Na I D, is strikingly different from the ionized gas (Figure 
[4|. At the 4(7 level, Na I D absorption extends from radii 
of 0.5 to 3 kpc, which is farther than the ionized outflow. 
(Closer than 0.5 kpc, the nuclear emission washes out 
the signature of the host galaxy.) A significant area of 
absorption is also seen atop the blue continuum peak 3.5 
kpc south of the nucleus. 

The neutral gas velocities are also much higher than 
the ionized gas velocities. The Na I D velocity maps 
consist of three distinct regions, whose velocity averages 
are found in Table [l] The lowest velocities lie atop the 
southern star-forming arc. Higher velocity, broader com- 
ponents are found in the outflow surrounding the nucleus 
(not including the northern quadrant). The highest ve- 
locities are found in the northern quadrant of the nu- 
clear outflow. In those parts of the northern quadrant 
that line up with the radio jet in Mrk 231, the velocities 
reach —1400 km s~^. 

The velocities observed in the arc are comparable to 
th ose measured in starburst-domina ted ULIRGs (Table 



l| ! Martin|[2QQ5| [Rupke et aL||2005c 

the nuclear wind are signiticantly h 



However, those in 

the nuclear wind are signiticantly higher; t hey are more 
like those found in som e Seyfert ULIRGs (Rupke et al. 
2005a||Krug et al.||2010[ ). 



4. DISCUSSION 
ULIRGs host powerful starburst-driven outflows 



( Heckman et al ||1990[ |2000| ) that arise in ah UL IRGs 
] Rupke et al.||2005c| aimsnow kpc scale extents (Mar- 
tin|| LiO06| |Shih fc Rupke] [2QTq| . They have {vqs%)~= 



—450 km s ^ in systems w hose infrared lu minosit y is 



dominated by star formation ( Rupke et al.|2005c^ ; Martin| 



2QQ5D . |Rupke et aP (| 2QQ5cD calculate that ULIRG mass 
outflow rates are 20% of the star formation rate on av- 
erage. However, evidence has remained elusive for large- 
scale out flows that are cl early AGN-driven in LIRGs or 
ULIRGs (Rupke et al. 2005a). 

The data we present here are clear: there is a neutral, 
~1000 km s~^ outflow in Mrk 231 that extends in every 
direction from the nucleus (as projected into the plane of 
the sky) out to at least 3 kpc. Such high velocities have 
not been seen in starburst ULIRGs (Table |T), providing 
strong circumstantial evidence that this wide-angle nu- 
clear wind is driven by radiation or mechanical energy 
from the AGN. 

Fro m what we know abo ut the structure of galactic 
winds fVeilleux et al. 2005), the molecular disk in Mrk 
231 collimates this nuclear wind. Given the disk's almost 
face-on orientation, we must be looking "down the bar- 
rel" of a biconical outflow. The other end of this bicone 
is receding from us, behind the galaxy disk and therefore 
invisible at optical wavelengths. 

It is apparent from the vel ocity map that the 
north-south r adio jet in Mrk 231 ICarilh et al |1998[ [uT 



vestad et al.|19 99 ) is coupling to the nuclear wind, accel- 



erating the neutral gas to even higher velocities. The jet 
is not constrained to emerge perpendicular to the disk, 
and thus produces an asymmetric effect. The present 
data imply that the northern arm of the large-scale jet is 
on the near side of the molecular disk. Neutral outflows 
driven by jet interactions with the interstellar medium 
(ISM) on kpc scales have also been observed in radi o 
galaxies in H I absorption ( Morganti et al.|[2005 , 2007). 
The newly discovered jet-wind interaction in Mrk 231 ap- 
pears to be similar, though this time the jet accelerates 
a n already in s itu wind. 
"Carilhet al.](|T998) and Taylor et al.|(|1999D studied the 



dittuse radio continuum emission from Mrk 231, which is 
symmetric about the nucleus on scales of 100 pc to 1 kpc. 
They hypothesized that this emission is produced by in 
situ electron acceleration, but could not rule out that 
the AGN distributes these electrons through a wide-angle 
outflow. Our data are further evidence for the in situ 
interpretation, since we now know that the AGN outflow 
reaches larger scales and is asymmetrically accelerated. 

Along with these AGN-driven outflows, Mrk 231 also 
hosts a starburst-driven wind. Blueshifted velocities with 



(^98%) 



-570 km s are observed in the star-forming 



arc south of the nucleus. These are comparable to the 
maximum velocities observed in other starburst ULIRGs 
(Table [T]). We cannot rule out that the underlying con- 
tinuum could be backlighting the nuclear wind at this 
location. However, a simpler explanation is that in situ 
star formation is driving an off-nuclear wind. 

The ionized outflow in Mrk 231 overlaps spatially with 
the neutral outflow near the nucleus and shows broad 
velocity proflles. It also shows higher velocities to the 
north. However, the velocities and extent of this wind 
phase are more modest. The ionized outflow in Mrk 231 
could be another phase of the nuclear wind, providing 
further evidence that the differ ent ^as phases in ULIRG 



winds are not s trongly coupled ( Rupke et al.|2005a Shih 
Rupke||2010| ). It may also be the classic AGN narrow 



line region (e.g., Crenshaw fc Kraemer|[2QQ5 ). 

The nuclear wind is clearly the on e detected previ ously 
in the neutral (Rupke et al. 2005C( ) and molecular (Fis- 



cher et al. 2010 Feruglio et al. 2010[) phases. Rupke 



et al. "i 2005cp found an extended wind ( on kpc scales) 



but on ly had data in a north-south slice. [Feruglio et al.| 
(2010) marginally resolved the molecular component out 
to radii of ~0.6 kpc, and estimated molecular mass and 
energy fluxes of several hundred Mq yr~^ and ~10^^ erg 

The current data are a significant improvement on 
these results. They clearly resolve the structure of the 
wind on sub-kpc scales and reveal that more than one 
power source is acting. To produce a wide-angle outflow, 
the AGN energy and momentum must be injected into 
the ISM over a substantial volume through mechanical or 
radiative processes. Either (1) wide-angle winds from the 
AGN accretion disk or broad line region or (2) coupling 
of the AGN radiation t o the dust in the wind are pos- 



sible explanations (e.g., Crenshaw et al. 2003 Hopkins 
^fc Elv is 2010). However, the high-velocity component of 
the wind that lies north of the nucleus must experience 
an additional, asymmetric force. The north-south jet in 
Mrk 231 is a natural explanation. 
Using a simple model, we can estimate the mass and 
fluxes (M and E) in the nuclear wind (Rupke] 



energy nuxes ( M and J m tne nuclear wma (nupKej 
et al ||2005c| |Shih fc Rupke||2010[ ). We assume a face-on 
wind for simplicity. We model the wind as emanating 
from the nucleus in a thin shell, and co mpute a time- 



averaged mass flux. Rupke et al. (2005c) also discuss a 



thick-shell model. However, computing the mass outflow 
rate for this model would require integrating hypothetical 
velocity proflles over each line of sight; such an exercise 
is beyond the scope of this Letter and would increase the 
number of unconstrained parameters. 

There are two signiflcant unknowns in the thin shell 
model. The flrst is the radius of the shell. Given the 
spatial smoothness of the Na I D velocity map, it must 
be larger than 2 kpc, or we would observe the veloc- 
ity to decrease (in projection) with increasing radius. 
However, it cannot be too large, or we would observe 
a more extended absorption region. We use a flducial 
radius of 3 kpc; larger values will increase the resulting 
M and E. The second uncertain ty is the ionization s tate 
of Na. As in previous work (e.g., Rupke et al.|2005b ), we 
use Milky Way measurements as a baseline and assume 
A/'(Na l)/A/'(Na) = 1 — y = 0.1, where y is the ionization 
fraction. The fluxes are inversely proportional to (1 — 7/), 
such that a more ionized wind will yield higher M and 

With this model, we compute M = 420 {R/3 kpc) 
[0.1/(1 - y)] Mq yr-i and ^ = 7.3 x 10"^^ {R/3 kpc) 
[0.1/(1 — y)] erg s~^. For comparison, the star for- 
mation rate and AGN luminosity in Mrk 231 are 172 
M0 yr~^ and 1.1 x 10^^ erg s~^. (These numbers are 
computed by assuming the bulk of the galaxy's radiation 
emerges in the far-infrared, using the AG N contribution 
to M rk 231 's bolometric luminosity from [Veilleux et aT] 
2009[ and cal culating the star formation rate with the 
formula from Kennicutt 1998 ) The wind is clearly re- 
moving signiflcant amounts oi gas from the nucleus - at 
a level of 2.5 times the star formation rate. This mass 
ou tflow rate is al s o stri kingly similar to that estimated 



b y [Feruglio et al. (2010) from a different tracer. 

I'he energy outflow rate is about 0.7% of the radia- 
tive luminosity of the AGN, meaning that only a small 
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part of the radiative output of the AGN has to couple to 
the outflow. This is remarkably similar to the required 
coupling efficie ncy needed for an AGN wind in the two- 



( 2Q1Q| ). In this model, 
, driving a diffuse out- 



phase model of Hopkins & Elvis 
the AGN couples to the hot iSlv:, 
flow. This diffuse outflow disrupts cold clouds, which 
are propelled outward by the diffuse wind and by radia- 
tion pressure from the AGN. Together, the hot and cold 
winds provide negative feedback on the AGN and star 
formation. 

5. SUMMARY 

The galaxy Mrk 231 is in the late stages of a major 
merger, hosts a QSO and an obscured starburst, and ex- 
hibits gas outflows in multiple forms. We show in this 
Letter that there is a massive, energetic, and neutral 
outflow extending in all directions from the nucleus (as 
projected in the plane of the sky) to radii of at least 3 
kpc. This outflow has blueshifted velocities that reach 
1100 km s~^. On the axis of the radio jet in Mrk 231, 
^98% is as large as 1400 km s~^. Mass and energy out- 
flow rates are ^2.5 x the star formation rate and >0.7% 
of the AGN luminosity. The former suggests strong neg- 
ative feedback to star formation, and the latter is con- 
sistent with the c oupling efficiency require d in the AGN 



feedback model of Hopkins & Elvis (2010) 



By comparing the velocities in this wind to those in the 
average starburst-driven ULIRG, we conclude that this 
wide-angle wind is likely powered by the central AGN, 
and its power source must be directed fairly symmetri- 
cally. The northern quadrant of this nuclear outflow ex- 
periences additional acceleration due to the radio jet in 
Mrk 231. Furthermore, a starburst-powered wind arises 
in the extended galaxy disk to the south of the nucleus. 

These data illuminate the structure and power source 
of the kpc scale wind in Mrk 231, and show that strong 
feedback occurs as a result. The ubiquity of winds 
in major mergers is well-known, but our understand- 
ing of their full impact on galaxy e volution is im proving 
through integral field observations ( Shih fc Rupke 2010 
and the present work), H erschel studies (Tischer et al. 



2010 Sturm et al. 



observations (ll^eruglio "et al.||2010[ ). 



2010 ), and gr ound-based molecular 
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Figure 1. Left: absorption line fits to Na I D. One spectrum is shown in each major outflow region (Sections |3.2| and [4]). The spaxel 
coordinates of each flt are shown in the upper left of each plot; the nucleus is located at spaxel [11,13], and spaxels are 0?^ in size. The 
vertical dashed lines show the locations of He I 5876 A and Na I D according to the rotation curve. The red line is the total (absorption 
+ emission) flt, while the blue line is absorption only. Right: emission line flts to the Ha / [N II] spectral region. The red line is the total 
(extended line emission + broad line region emission + continuum) flt, while the blue line is broad line region emission + continuum only. 



Table 1 

Outflow Velocities in Mrk 231 



Label 


Phase 


(FWHM) 


Ko%> 


(^84% > 


(^98% > 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Mrk 231 starburst wind 


Neutral 


510 


-150 


-360 


-570 


Mrk 231 nuclear wind 


Neutral 


600 


-360 


-620 


-880 




Ionized^ 


690 


-160 


-450 


-750 


Mrk 231 jet-influenced region 


Neutral 


640 


-520 


-800 


-1070 


Starburst ULIRGs^ 


Neutral 


330 


-170 


-310 


-450 



Note. — Column 2: gas phase. Column 3: full width at half maximum of 
velocity distribution, in km s"-*^. Column 4: velocity at the center of (Gaussian) 
velocity distribution. Column 5: '^^84% = '^50% ~ §4% of the outflowing gas has 
velocities closer to 0. Column 6: ^'98% = '^50% ~ 2cr. 
^ For the broad, outflowing component. 

^ Computed from starburst-driven ULIRGs in Rupke et al.| ( 20Q5c| . This is a 
sample average, as opposed to the spatial average tor Mrk 2'6l. 
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Figure 2. Top: continuum image of Mrk 231 in the 435W and 814W filters, from Hubble Space Telescope observations with the Advanced 
Camera for Surveys. The field of view of our CMOS data (6(^3 x 7^5) is overlaid as a box. Bottom left: host galaxy continuum image in 
our CMOS field, summed over the wavelength range 6100 — 6900 A, in logarithmic fiux units. Bottom center: Ho; emission, in logarithmic 
flux units. Bottom right: map of the [N II]A6583/Hq; flux ratio. 
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Figure 3. Left: central velocity and FWHM maps of the rotating ionized gas component, in km s"-*^. Th e ionized gas rotation is 
consistent with that of the nuclear molecular disk. The line of nodes position angle is from Downes Sz Solomon] |l998 ). Right: ^'50% and 
VQg% (Section[2| for the broad, outflowing ionized gas component. A near-nuclear outflow is present, as well as a highly blueshifted H II 
region in the southwest. 
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Figure 4. Equivalent width, central velocity, FWHM, and '^^98% maps of Na I D. A nuclear outflow extends from the nucleus up to 
2 — 3 kpc in all directions (as projected in the plane of the sky). The high velocities suggest that the AGN powers the nuclear wind. The 
northern quadrant of the nuclear wind is further accelerated by the radio jet. A lower- velocity starburst-driven outflow is present in the 
south. 



